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Abstract 

We report here the existence of Ermanno-BernouUi type invariants 
for the Manev model dynamics which may be viewed upon as remnants 
of Laplace- Runge-Lenz vector whose conservation is characteristic of the 
Kepler model. If the orbits are bounded these invariants exist only when 
a certain rationality condition is met and thus we have superintegrability 
only on a subset of initial values. We analyze real form dynamics of the 
Manev model and derive that it is always superintegrable. We also discuss 
the symmetry algebras of the Manev model and its real Hamiltonian form. 

1 Introduction 

Since time immemorial the circular motion was the archetype motion of the 
heavenly bodies, and circle was assumed to be embodiment of perfection — in the 
East and the West alike. No wonder that when the observation data challenged 
this view a superposition of several circular motions or circular motion with an 
off-centre Sun were proposed in order to keep the circular paradigm intact. 

Since Kepler and Newton elliptical trajectories replaced circular ones as an 
archetype of the (bounded) planetary motion and the circle is nowadays viewed 
upon rather as a degenerate ellipse than as an ideal incarnated. The archetype 
of elliptical motion is even exported to the atomic realm as we see depictions of 
'the atom' with ellipses representing electrons' motion around the nuclei. 

The advent of Einstein's theory did not produce a new archetype of heavenly 
motions, apart from the exceptional case of a collapse into the (still hypothetical) 
black holes. Nevertheless, among the variety of relativistic effects the perihelion 
advance of inner planets is definitely the best recognizable effect in the Solar 
system. Maybe it is time to accept a new archetype: precessing ellipse (or 
more generally, precessing conies). There are also classical arguments in its 
favour: Kepler-type motion is generally not preserved by small perturbations 
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and generally any sort of 'real world' interactions like Solar pressure, drag etc 
would destroy 'fixed ellipse' motion If precessing conies give us 'the typical' 
motion of planets it is tempting to ask which central force field produces them. 
Surprisingly or not, the answer is: the Manev model j2|- Here we already have 
persistent KAM tori and cylinders for a large class of even non Hamiltonian 
perturbations |T] and this is an additional argument in favour of the Manev 
model. 

Kepler problem is famous as one of archetypes of superintegrable systems 
and probably, the first one where an unexpected non-Nother symmetry has 
been uncovered. It is intriguing to ask whether Manev problem shares this 
property and here we report that this is indeed the case, but not for all initial 
data. Artificial models presenting such a behaviour are already known, e.g. 
sum of squares of Hamiltonians of independent oscillators, but Manev model 
is a 'real world' example having this property. Let's remark that for a generic 
central potential we could have disjoint set of initial data corresponding to closed 
orbits but in our case all points on certain level sets of the angular momentum 
lie on closed orbits which are intersections with the level sets of the additional 
invariant. 

Also we will show that real form dynamics of the Manev problem — a closely 
connected dynamical model which we shall introduce below — is superintegrable 
for all initial data. Real form dynamics of the Manev problem is interesting 
enough to validate a separate analysis and we will describe it at the end of the 
article. 



2 The Manev Problem Basics 

By Manev model |3] we mean here the dynamics given by the Hamiltonian: 

H=lipl+Pl+Pl)-^-^ (1) 

where r = ^yx^ + y^ + z^ ■, A and B are assumed to be arbitrary real con- 
stants whose positive values correspond to attractive forces. The genuine model 
proposed by George Manev was not invented as an approximation of relativ- 
ity theory but as a consequence of Max Planck's (more general) action-reaction 
principle and it derived a specific value for the constant B = Nevertheless, 
Manev model offers a surprisingly good practical approximation to Einstein's 
relativistic dynamics — at least at a solar system level — capable to describe both 
the perihelion advance of the inner planets and the Moon's perigee motion. In 
the last decade it had enjoyed an increased interest either as a very suitable 
approximation from astronomers' point of view or as a toy model for applying 
different techniques of the modern mechanics (see e.g. 01 El 13 El ) ■ 

Due to Hamiltonian's rotational invariance each component of the angular 
momentum 

Lj^e-jkmPkXm {Lj ,Lk} = e-jkni Lrn with {xi,X2,X3) ^ {x,y,z) (2) 
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is an obvious first integral: {H, Lj} — and so, like the Kepler problem, the 
Manev model is integrable. Components themselves are not in involution but 
span an so(3) algebra with respect to the Poisson brackets. 

The dynamics is confined on a plane which we assume to be Oxy and is 
separable in radial coordinates. On the reduced phase space (see e.g. [5| for the 
generalities of the reduction procedure) obtained by fixing the angular momen- 
tum = L to a certain value £ the motion is governed by: 

The dynamics behave like radial motion of Kepler dynamics with angular mo- 
mentum squared £^ — 2B; while the case 2B > £^ corresponds to overall cen- 
tripetal effect. On the other hand, the angular equation of motion 6 — L/r^ 
is still governed by the 'authentic' angular momentum i (and r is as just de- 
scribed). Consequently, the remarkable properties of Kepler dynamics that all 
negative energy orbits are closed and the frequencies of radial and angular mo- 
tions coincide (for any initial conditions) are no more true. Thus we may have 
not only purely classical perihelion shifts but also if 2B > £^ ^ we may have 
collapsing trajectories which are spirals, even though in phase space they are 
symplectic transformations; while in the Kepler dynamics the only allowed fall 
down is along straight lines. For this reason the set of initial data leading to col- 
lision has a positive measure and this may offer an explanation why collisions 
in the solar system are estimated to happen more often than Newton theory 
predicts [TU] . 

The dynamics of Manev model has already been thoroughly analyzed (see 
e.g. El El 13 E]) and we shall concentrate here on some of its invariants and 
resulting symmetry algebra. 



3 The Kepler Problem Invariants 

In the case of Kepler problem, corresponding to i3 = 0, we have more first 
integrals (for details and historical notes see e.g. [TTl IT^ IT^ ITlj l: 

A A ^ 

J^=PyL-\ X, Jy^-pxL^ y, {Hk,J} = 0. (4) 

r r 

where Hk is the Kepler Hamiltonian and and Jy are the components of the 
Laplace-Runge-Lenz vector. They are not independent since 

J"^ = 2HkL^ + A^. (5) 

Together with the Hamiltonian and angular momentum they close on an algebra 
with respect to the Poisson brackets: 

{Hk,L}^Q, {i,Jj = Jj,, {L,Jy} = -J,, {J^,Jy} = -2HKL. (6) 
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After redefining E = J j \/\ — 2Hk\ we get: 



{L, Ex} = Ey 



{L, Ey] — —Ex 



{Ex,Ey}^-sign{HK)L (7) 



which makes obvious the fact that we have an so(3) algebra for negative energies 
and so(2, 1) for positive ones. In the case of the 3-dimensional Kepler problem 
the components of the angular momentum give us another copy of so(3), see 
eq. l(2Jl, so the full symmetry algebra is so(4) or so(3, 1) depending on the sign 

According to the first use of these first integrals was made by J. Hermann 
(= J. Ermanno) in 1710 (in order to find all possible orbits under an inverse 
square law force) in the disguise of 'Ermanno-BernouUi' constants: 



J i: — Jx ~F l"-^y — 



— A ^ i Lpr 



which satisfy: 

{Hk,J±} 



0, 



{L,J±} = ±iJ±, {J+,J-} 



-AiHkL. 



(8) 



(9) 



Curiously enough, the initial Kepler idea for a circular orbits with an off- 
centre Sun happens to be correct in a different context. As discovered by Hamil- 
ton '16', the velocity vector in the Kepler dynamics moves along a circle laying 
in a plane containing the origin but, in general, not centered, at the origin. If 
we choose the a;-axis pointing to the point of closest approach, the centre of this 
circle is located at (0, eA/£, 0) and its radius is A/i where e = -^71 + 2EP / A^ is 
the eccentricity of the orbit. If the spatial orbit is hyperbola (e > 1) the velocity 
space orbit is only a section of the circle, otherwise the full circle is traversed. 

4 The Manev Problem Invariants 



Assuming Q ^ I > 2B and denoting 
^ f- - 2B 



w 



L2 

r 



one easily verifies that 



_d 



vw ± i—w I e 
dt^ 



±ivB 



. 2 . 

and since = we obtain: 

d6 L at 



d_ 

di 



P - 2B 



A 



= 



0. 



(10) 



(11) 



(12) 
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In the case when ^ 7^ 0, > 2B , H < and A > the motion is on a 
2-dimensional torus. In order to have globally defined constants of motion in 
this case we have to require that the ;/'s be rational i.e. 



\/p - 2B : e = 



k 



(13) 
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with m and k integers. Then due to eq. 

m k 



J±=J. 



-A =F iipr 
k r m 



are conserved by the flow of eq. iQ]) on a surface L — £ satisfying the rationahty 
condition Thus we have conditional constants of motion which exist only 

for disjoint but infinite set of values £ (c.f. invariant relations of |15j). 

The trajectory in the configuration space is a 'rosette' with m petals and 
this is connected to the fact that J7± are invariant under the action of the cyclic 
group generated by rotations by angle 2^ : 

k 

e-^9 + 2TT—n n = 0,l,..., 771-1. (15) 

777 

While in the Kepler case we could unambiguously attach the Laplace-Runge- 
Lenz vector to Ermanno-BernouUi invariants this is not possible now due to this 
finite symmetry. (It is intuitively clear that if the Laplace-Runge-Lenz vector 
points to the perihelion of the Kepler ellipse, now we have m petals to choose 
between.) Anyway, up to this unambiguity, or restricting ourselves to one of 
the 777. sectors we may note that while the radial/angular components of the 
Laplace-Runge-Lenz vector take the form: 

^2 

Jr = A Je = -Lpr (16) 

r 

in our case Jr + iJe = ( J^— — — — iipr) e'^"^"^^^ and hence: 



£^ A" 

Jr = [v ) cos(t/- 1)61 + ^p^sin(7/- 1)61 (17) 



r V 



P A\ 

Je — —£pr cos(iy — 1)6* + ( 7/ I sin(i^ — 1) 



Turning to the algebraic properties of the new invariants one finds that the 
Poisson brackets between the real and imaginary parts of J± = Jo T ^Ji = 
yJx "F iiTy are: 



{if, Jo,i} 



0, 



m 



{i,Ji} = -^Jo 



{Jo,Ji} = -^L 



2H 



2B f2L^ 



(18) 



Here we have a 1/r^ term which seems to obstruct the Poisson brackets to 
form a closed algebra. Fortunately, redefining 5o,i = Jo,il L^/TF^^l? we obtain 
the closed algebra Qh.l'- 

TD 777 

k k 
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k (L2-^2)2 



2HL 



e A^ 2L^ - £^ 



(19) 
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in which iJ is a central element and L, £q and Si can be viewed as Cartan and 
root- vector generators. Due to l(T^ Qh.l is a deformation of gl{2). Of course, 
we have in addition the so(3) algebra l(2Jl- 

Similarly, in the case when Q ^ (?■ < 2B we may denote '^'^^^ — with v 
real and 



A 



(20) 



are first integrals for any I and they satisfy: 



{i/,f±} = 0, 



{L,£±} = TvS± 



2v 



(L2 _ £2)2 



2HL 



A^2I?_f_ 



(21) 



The algebra g'^ ^ satisfied by i7, L and £± is quite analogous to Qh.l but with 
a different function at the right hand side of the bracket {£+,£-}. 
Finally, when t"^ = 2B we have the first integral: 



j = Lpr + AO 



(22) 



satisfying {H,j} ^ 0, {i, j} = A. 



5 Real Form Dynamics 

Here we briefiy recall the notion of real form (RF) dynamics referring the reader 
to J7j for more details and a list of examples. 

We start from a standard (real) Hamiltonian system 7i = {A^,!^,//} with 
n degrees of freedom and at the present stage we assume that our phase space 
is just a vector space M — M^", 

Let's consider its complexification: US' = { Al^, i/'', cj*'} where M!^ can be 
viewed as a linear space over the field of complex numbers: 

In other words the dynamical variables in M!^ now take complex values. We 
assume that the Hamiltonian H (as well as all other possible first integrals in 
involution Ik) are real analytic functions on M. which can naturally be extended 
to M!^. We introduce on the phase space M an involutive, symplectic auto- 
morphism C : M ^ M: 

C^^l, C{{F,G}) = {C{F),C{G)} (23) 

where with some abuse of terminology we use the same notation for the action 
of C on the dual of the phase space. 

Since C has eigenvalues 1 and —1, it naturally splits A4 into two eigenspaces: 

M^M+®M^ (24) 
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whose dimensions need not be equal. Due to the fact that C is symplectic 
and M+ are symplectic subspaces of M and we will write = cj+ © 
Assuming a symplectic frame adapted to C we have: 

w = ^ Apk+ A d(7fe+ + ^ dpfe_ A d(7fe_ . 
fc=i fc=i 

The automorphism C can naturally be extended to M."^ and it splits it again 
into a direct sum of two eigenspaces: 

Similarly, the action of the complex conjugation * produces splitting into real 
and imaginary parts of the corresponding spaces. By construction C commutes 
with * and their composition C = C o * — * o C is also an involutivc symplectic 
automorphism on then we define AIr to be the fixed point set of C i.e. 

and it is again a symplectic subspace. From now on wc will be interested in 
dynamics on and its connection to the initial real dynamical system. 

In order to construct 'real form dynamics' we shall assume that the Hamil- 
tonian is C-invariant, i.e.: 

C[H) ^ H. (25) 

Then the Hamiltonian on the complexified phase space (being the same 
analytical function of the complexified variables) will share this property. 
The 'real form dynamics ' may be defined either as: 

i) complexified Hamilton equations on Ai"" being consistently restricted to AIr. 

This gives a real vector field tangent to Mm and satisfying the equations 
of motion given by the real part of or 

ii) dynamics on Mr defined by the restricted H'^ and w'' (whose restrictions 

are real on Mr): 

^1 ^ H + C(H)^H + CiHr 



(jj 



I Mm 



^ dRep^ji AdReg'j: -dlmp^ Adlmg^. (26) 



Now we have a well defined dynamical system Hm — {Mm, i^Itv^j,, H\j^^} 
with real Hamiltonian and real symplectic form on a subspace of the complexified 
phase space. 

It is noteworthy that the 'real form dynamics' corresponding to a Liouville 
integrable Hamiltonian system is Liouville integrable again 17,. Similarly, the 
'real form dynamics' corresponding to a superintegrable Hamiltonian system is 
superintegrable again. In such a case we have n € [n + 1, 2n — 1] independent 
constants of motion which are no more in involution. It could easily be checked 
that they will again produce k independent constants of motion of the RF 
dynamics. 
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6 Real Form Dynamics of the Manev Problem 



The Manev Hamiltonian (and the canonical symplectic form as well) is invariant 
under the involution C reflecting the y-degree of freedom: 

C{x)=x, C{y)^-y, C{z) = z 

C{Px)^Px, C{py) = -py, C{pz)=Pz- (27) 

Consequently, the 'real form dynamics' of Manev model for this choice of invo- 
lution will be given by: 

Hr = lipl-pl+pl)---^ (28) 
= dpx A dx — dpy A dy + dpz A dz 

where p = i/a;^ — + is the 'radius' of the pseudo-sphere . This is not an 
ordinary central field dynamics but rather an 'indefinite metric central field' as 
-Hr depends on indefinite metric distance p. The real form Hamiltonian H-r 
and the appropriate 'angular momentum' Lj are still commuting first integrals 
and the model is integrable. The involution acts on Lj according to: C{Lj) = 
{-lyLj and 

{^■,Lfc} = e,fe,(-l)^+'^+iZ, (29) 

instead of eq. |5J); the corresponding algebra is so{2, 1) which is the real form of 
so(3) obtained with a C-induced Cartan involution. 

We shall assume again that the motion is on the Oxy-plane and in order to 
avoid the question of the behavior of trajectories on the singularities we restrict 
our attention on the C-invariant configuration space: 

{{x, y,z) eM^ \ z = 0,x^ > y^, x > 0}. 

Then the dynamics is separable in pseudo-radial coordinates = artanh(y/a:) G 
(— oo, oo) and p G (0, oo): 



UJ = dpp Adp + dvr^ A di? 

with L = Lz ~ TTij, hence tt'^ = and i9 = —L/p^. Due to the different 
symplectic form L generates now transformations which preserve p. 

The type of the p-trajectories could be easily read off after the observation 
that the value of the real form Hamiltonian will be: 

1 / .2 A B 



2 V P J P P 

due to x^ — y^ = fP' — j p^ and denoting the value of L by I. Introducing v — p 
and u — \/ p we obtain an equation describing conies in the {u, u)-space: 

u^{l^ + 2B) + 2Au + {2h - v^) = 0. 
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Performing the same type of analysis as in p] we may conclude that we may 
have three types of qualitatively different dynamical regimes: 



for l"^ + 2B > we have a family of hyperbolas. 

for P + 2B = we have a family of parabolas for A ^ which degenerate at 
yl = into pair of lines parallel to the l/p-axis. 

for £^ + 2B < (only possible for repulsive Manev term) we have a family of 
ellipses. 

Of course, in all these cases we have to exclude the region u < 0. 

In order to obtain more specific information about the motion we will need an 
equation for the trajectories. Let's note again that in the case of non-vanishing 

jO-motion takes the form: 



angular momentum we have dt = —Ij-dd. As a result the equation for the 



di?2 p p 
Assuming + ^ we introduce 

e + 2B P 



d2 L2 ^2^25^2 

A = 0. (31) 



and obtain an inverted oscillator-type equation: 

d2 



w-v'^w = 0. (33) 



Denoting by Cj the integration constants below we conclude that: 
If ^2 _^ > the solution p"^(z9) will be: 

Z2 A 

— = ci cosh(wi9) + C2 sinh(i;'!9) ^. (34) 

p 

Trajectories may collapse (p — + 0) for d — + ±oo and ci > C2 > 0, or p may 
tend to oo as z9 tends to certain values 'dmin and i)max- 

lil"^ + 2B^Q the solution of eq. (EB wiU be: 

p-^ = ^ ^ (35) 

If A > trajectories collapse for -d — > ±oo and if A < then p ^ oo as 
■d tends to some t?mm and d„iax- The case A — Q leads to linear solution 
p~^{'d) and corresponds either to the motion along fixed p or to trajectory 
starting at p = cxd and some value of and collapsing for § oo (or its 
reverse) . 
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If. 



^ + 2B < we will have a solution which oscillates harmonically between 
some values prmn and Pmax- 



12 ^ 

— = C5 coslvt)) + cq sinivd) tt- 

p v' 



(36) 



The case when pmin < means that 'acceptable' motions will be trajec- 
tories coming from p = 00 at d = "dmin and going to p = cxd for some 

In the special case of vanishing angular momentum we have 1-dimensional 
motion along the ray d = const. It may be oscillating between some Pmin and 
Pmax or heading to collapse, or to infinity. 

It is worth noting that the only trajectories which are compact in the {x, y)- 
space are those collapsing at both their ends in the origin tangentially to the 
boundaries x — ±y and those oscillating on a line interval of -d = const. This 
is to be contrasted to the standard Manev or Kepler problems. 

Obviously, when B = we will obtain a real form dynamics of the Kepler 
model. In this case we have even fewer possibilities for compact trajectories as 
we could not have oscillations along the line of i? = const. 

Since the motion is never on a 2-torus the additional first integrals are always 
globally defined for all initial data. When 7^ £^ > —2B they take the form: 



J± = 



A 



± Lpp 



(37) 



As we did earlier we can introduce the renormalized E± = J±IL\J — P and 
derive for them the following symmetry algebra q'^ ^ 



{ffR,f±} = 



2v 



(i2 _ 



2HrL- 



2L''-e 



(38) 



Like in (|19|l above q'^ ^ deformation of gl{2) having the same i7, L depen- 
dence in the right hand side of 1)38(1 . though L and L have different properties. 

Note that the algebras Qh,l and q'^ ^ seem very close, i.e. they do not 
change effectively when passing from one real Hamiltonian form to the other. 
The reason for this is the fact, that all its generators are invariant with respect to 
the involution C. The situation changes when we consider the algebra satisfied 
by Lj , see eq. (j^HJ. 



In the case when Q ^ i"^ < -2B let v"^ 



gi^^^ and then invariants are: 



J± 



Jo T i Ji 



i2 A ~ 

V =F iLpp 

p V 



(39) 
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Redefining again: £0.1 = Jq^\IL\/lP^--P- wc obtain the brackets: 



{iJK,fo,i} — 0, 



(40) 



V 



(L2 _ py 



V 



Finally, when = 2B we have the first integral: 



j = Lpp - Ad 



(41) 



satisfying {H,j} = 0, {L,j} = -A. 

7 Conclusions 

The discovered existence of Ermanno-BernouUi type invariants strengthens our 
belief that Manev model has an exceptional position among the central field the- 
ories. Not only it provides a better description of the real motion of the heavenly 
bodies than Kepler model but to a large extent it shares its superintegrability — 
probably its most celebrated mathematical feature. As a result it provides 
also a testbed for analysing the intricate interplay between integrability and 
superintegrability — a testbed having the advantage of being realistic and intu- 
itive interaction. 

Also, from the viewpoint of a RF dynamics enthusiasts we see here a curi- 
ous (and encouraging) example when the RF dynamics — exotic as it may be — 
behaves 'better' than the original problem as it is always superintegrable. 

As a final remark we note that the new results reported here do not contra- 
dict existing classifications of superintegrable models (see e.g. ^Hl) where only 
first integrals which are low order polynomials in momenta are described. 
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